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We  present  atomistic  first  principles  results  for  the  lattice  thermal  conductivity  of  GaN  and  compare 
them  to  those  for  GaP,  GaAs,  and  GaSb.  In  GaN  we  find  a  large  increase  to  the  thermal  conductivity  with 
isotopic  enrichment,  ~65%  at  room  temperature.  We  show  that  both  the  high  thermal  conductivity  and  its 
enhancement  with  isotopic  enrichment  in  GaN  arise  from  the  weak  coupling  of  heat-carrying  acoustic 
phonons  with  optic  phonons.  This  weak  scattering  results  from  stiff  atomic  bonds  and  the  large  Ga  to  N 
mass  ratio,  which  give  phonons  high  frequencies  and  also  a  pronounced  energy  gap  between  acoustic  and 
optic  phonons  compared  to  other  materials.  Rigorous  understanding  of  these  features  in  GaN  gives 
important  insights  into  the  interplay  between  intrinsic  phonon-phonon  scattering  and  isotopic  scattering  in 
a  range  of  materials. 
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Introduction. — Gallium  nitride  (GaN)  is  a  wide  band  gap 
semiconductor  and  a  promising  candidate  for  use  in  opto¬ 
electronic  devices  [1,2],  high-frequency  switches  [3,4], 
and  high-power  electronics  [5].  Thermal  management  is 
of  critical  importance  in  such  devices,  especially  as  typical 
device  sizes  continue  to  decrease.  GaN  has  the  attractive 
feature  of  an  unusually  high  lattice  thermal  conductivity, 
kl,  with  measured  room  temperature  values  around 
230  Wm_1  K_I  [6-8].  To  date,  the  origin  of  its  high  kl 
is  not  well  understood  at  a  quantitative  level.  Furthermore, 
increase  of  its  kl  could  lead  to  better  performing  and 
longer-lived  electronic  devices. 

As  in  other  semiconductors,  heat  in  GaN  is  carried  pri¬ 
marily  by  phonons.  The  kl  is  limited  by  intrinsic  phonon- 
phonon  scattering  around  room  temperature  and  higher  [9] 
and  by  extrinsic  scattering  due  to  point  defects  such  as 
isotopes.  Recent  work  has  estimated  that  the  room  tempera¬ 
ture  “isotope  effect”  given  by  P  =  100  X  (/cpure/  ^natural  —  1) 
is  a  relatively  small  4%-15%  in  GaN  [6,10-12].  Here  /cpure 
corresponds  to  kl  for  isotopically  pure  GaN  and  /cnatural  to 
kl  with  naturally  occurring  Ga  isotope  concentrations. 
Other  work  has  suggested  that  kl  of  defect-free  GaN  is 
much  larger  than  measured  values  but  that  the  enhancement 
comes  from  removal  of  vacancies  and  substitutional  defects 
rather  than  from  isotope  enrichment  [13,14]. 

Reliable  determination  of  k,  and  the  effects  of  different 
scattering  mechanisms  depends  on  an  accurate  representa¬ 
tion  of  the  anharmonic  phonon-phonon  scattering.  In  simple 
models  for  kl,  this  scattering  is  often  estimated  by  using  a 
mode  averaged  Gruneisen  constant  and  using  parameters 
adjusted  to  fit  measured  data  [10,11].  Such  models  lack 
predictive  capability.  Alternatively,  molecular  dynamics 
simulations  have  been  used  [15],  but  these  cannot  represent 
quantum  mechanical  scattering  and  so  are  of  uncertain 
validity  below  the  Debye  temperature  of  a  material 
(—600  K  for  GaN).  In  this  Letter,  we  employ  a  first 


principles  approach  to  definitively  determine  kl  in  GaN 
due  to  phonon-phonon  and  isotopic  impurity  scattering. 
Unlike  previous  estimates,  our  state-of-the-art  calculations 
show  that  the  isotope  effect  in  GaN  is  very  large,  —65%  at 
room  temperature,  making  it  comparable  to  the  highest 
isotope  effects  observed  in  diamond  (—50%)  [16],  graphene 
(60%)  [17],  and  boron  nitride  nanotubes  ( —  50%)  [18].  As  a 
result,  we  find  that  the  intrinsic  upper  limit  to  kl  in  GaN  is 
around  400  W  m~ 1  K  1  at  300  K,  considerably  higher  than 
currently  measured  kl  values.  We  find  that  these  features 
are  a  consequence  of  (i)  the  substantial  isotope  mixture  in 
Ga,  (ii)  the  large  energy  gap  between  the  acoustic  and  optic 
phonons,  and  (iii)  the  high  frequency  scale  of  the  phonon 
dispersion  in  GaN.  The  latter  two  points  weaken  the 
phonon-phonon  scattering  and,  in  particular,  the  scattering 
between  the  heat-carrying  acoustic  phonons  and  optic  pho¬ 
nons.  This  leads  to  higher  kl  and  makes  isotope  scattering 
relatively  more  important.  We  illustrate  this  effect  by  com¬ 
paring  the  kl  for  GaN  with  that  for  GaAs,  GaP,  and  GaSb. 

Ab  initio  thermal  transport. — We  employ  a  first  prin¬ 
ciples  approach  to  calculate  kl,  which  combines  an  exact 
numerical  solution  of  the  Peierls-Boltzmann  transport 
equation  with  accurate  calculations  of  the  harmonic  and 
anharmonic  interatomic  force  constants  (IFCs)  using  den¬ 
sity  functional  theory  and  density  functional  perturbation 
theory.  The  absence  of  accurate  values  for  the  anharmonic 
IFCs  has  limited  quantitative  results  for  kl  and  full  under¬ 
standing  of  the  underlying  physics.  A  similar  approach 
previously  used  by  us  has  accurately  described  thermal 
transport  in  silicon,  germanium,  and  diamond  [19,20]. 
We  consider  kl  limited  by  intrinsic  (three-phonon)  and 
isotope  defect  scattering  processes.  The  thermal  conduc¬ 
tivity  tensor  is 

«ap  =  X  J (dnV ^ \pd0-  (1) 
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In  Eq.  (1),  coA  and  vKa  =  dcoA/dqa  are  the  frequency  and 
ath  component  of  the  velocity  for  a  phonon  in  mode 
A  =  (q,  j),  respectively,  with  wave  vector  q  in  branch  j, 
nA  is  the  Bose  distribution  function,  and  tAq.  is  the  phonon 
lifetime.  GaN  has  a  wurtzite  structure  with  an  underlying 
hexagonal  lattice.  We  define  the  layers  to  be  in  the  x-y 
plane  with  the  c  axis  in  the  z  direction.  Then  Ka/3  is 
diagonal  with  one  out-of-plane  component,  kzz,  and 
two  in-plane  components,  kxx  and  Kyy.  In  principle,  for 
infinite,  perfect  hexagonal  systems  kxx  =  Kyy,  however, 


isotope  scattering  introduces  in-plane  anisotropy.  At 
around  room  temperature  and  higher,  this  anisotropy  is 
only  a  few  percent,  so  we  define  Km  =  (kxx  +  Kyy)/ 2  and 
/r0ut  =  KZz  to  discuss  in-plane  and  out-of-plane  thermal 
transport. 

The  linearized  Peierls-Boltzmann  transport  equation  is 
solved  by  using  an  iterative  scheme  described  in  detail 
previously  [19-21].  The  intrinsic  three-phonon  scattering 
probabilities  are 
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Here,  N0  is  the  number  of  unit  cells,  and  the  ±  correspond 
to  the  three-phonon  processes  that  satisfy  conservation  of 
energy  and  momentum:  coj(q)  ±  coji(q')  =  C0jii(q")  and 
q  ±  q'  =  q"  +  K,  where  K  is  a  reciprocal  lattice  vector, 
which  is  zero  for  normal  processes  and  nonzero  for  um- 
klapp  processes  [9].  Also,  Oai0r(OK,  I'k',  l" k")  are  third- 
order  anharmonic  IFCs,  extK  is  the  ath  component  of  the 
phonon  eigenvector  eA  for  the  Kth  atom  of  the  unit  cell  in 
mode  A,  and  Rl  is  the  lattice  vector  locating  the  Zth  unit 
cell.  The  scattering  time  due  to  isotopic  impurities,  t™,  is 
given  by  [22] 


1/t‘aS°  =  ‘  *** I25("a  “  "a')-  (3) 

*  A' 


In  Eq.  (3),  gK  =  ^ ~  mK)2  is  a  mass  variance 
parameter  with  /%  and  miK  being  the  concentration  and 
mass  of  the  ith  isotope  of  the  wth  atom,  respectively,  and 
mK  being  the  average  mass.  For  gallium  (60.11%  69Ga 
and  39.89%  71Ga),  gGa  =  1.97  X  10-4,  while  for  N,  the 
isotope  variation  is  negligible,  gN  =  0. 

The  formalism  described  above  requires  harmonic  IFCs 
to  calculate  phonon  frequencies  and  eigenvectors  from 
diagonalization  of  the  dynamical  matrix.  It  also  requires 
the  Oa;gr(0/c,  I'k1,  l" k")  to  allow  calculation  of  the  matrix 
elements  for  phonon-phonon  scattering  in  Eq.  (2).  These 
IFCs  were  calculated  within  the  framework  of  density 
functional  theory  and  density  functional  perturbation  the¬ 
ory  using  norm-conserving  pseudopotentials  in  the  local 
density  approximation  (LDA)  with  the  QUANTUM 
ESPRESSO  package  [23].  The  ground  state  configuration 
for  the  wurtzite  lattice  was  determined  by  varying  the  in¬ 
plane,  a ,  and  out-of-plane,  c,  lattice  constants  to  find  the 
minimum  energy,  and  the  internal  parameter  u  was  chosen 
so  that  the  residual  Hellman-Feynman  forces  on  the 
atoms  were  zero.  We  find  lattice  constants  that  are  1.8% 
less  than  experiment  (oexp  =  3. 190  A,  cexp  =  5. 189  A, 


and  wexp  =  0.377  [24]).  This  is  consistent  with  a  known 
shortcoming  of  LDA  approaches,  which  overbind  atoms  by 
l%-2%  [25].  Inclusion  of  zero-point  and  finite  temperature 
atomic  motion  [26]  into  LDA  results  gives  an  increase  of 
only  about  0.2%  to  the  zero  temperature  lattice  constants. 
Since  the  values  of  lattice  constants  affect  vibrational 
properties,  it  is  important  to  correct  for  this  limitation  of 
LDA  calculations.  Here  we  give  results  using  harmonic  and 
anharmonic  IFCs  calculated  for  the  following  two  choices: 
set  1  a  =  3.138  A,  c  =  5.109  A,  and  u  =  0.377  deter¬ 
mined  by  energy  minimization  with  a  0.2%  increase, 
and  set  2  a  =  3. 164  A,  c  =  5. 150  A,  and  u  =  0.377  deter¬ 
mined  by  energy  minimization  with  an  increase  of  1 .0%  in 
the  lattice  constants,  more  consistent  with  measured  values. 
Figure  1  shows  the  calculated  phonon  dispersions  for  wurtz¬ 
ite  GaN  using  set  1  (dashed  lines)  and  set  2  (solid  lines) 
compared  with  experimental  inelastic  x-ray  data  (circles)  in 
the  high-symmetry  directions  [27].  For  the  low-frequency 
acoustic  branches  there  is  little  difference  between  the 
two  sets.  However,  the  high-lying  optic  phonon  branches 
are  shifted  to  lower  frequencies  with  the  larger  lattice 
constants.  These  shifts  are  consistent  with  the  ab  initio 
calculations  in  Ref.  [27],  which  scaled  the  LDA  phonon 
dispersion  by  0.97  to  fit  the  experimental  data  better. 

To  determine  the  anharmonic  IFCs,  we  used  a  real-space 
approach  similar  to  that  of  Ref.  [28].  Pairs  of  atoms  were 
systematically  perturbed  from  the  ground  state  configura¬ 
tion  of  a  large  supercell.  The  resulting  Hellman-Feynman 
forces  on  all  the  atoms  were  calculated  via  T-point  self- 
consistent  calculations  for  a  number  of  different  perturba¬ 
tions  to  obtain  numerical  derivatives  of  the  forces,  which 
gave  the  anharmonic  IFCs  for  three-phonon  scattering.  For 
wurtzite  GaN  we  used  108  atom  supercells  and  calculated 
anharmonic  IFCs  involving  a  unit  cell  atom  and  its  27 
nearest  neighbors  (4th  neighbor  shell).  We  also  considered 
the  less  common  metastable  zinc  blende  phase  of  GaN  for 
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FIG.  1.  Calculated  phonon  dispersion  for  wurtzite  GaN  in 
high-symmetry  directions  with  set  1  (dashed  lines)  and  set  2 
(solid  lines)  lattice  constants  described  in  the  text.  The  hollow 
circles  correspond  to  experimental  data  determined  by  inelastic 
x-ray  scattering  [27]. 

which  we  used  216  atom  supercells  and  calculated  anhar- 
monic  IFCs  for  a  unit  cell  atom  and  its  29  nearest  neighbors 
(3rd  neighbor  shell).  In  principle,  anharmonic  terms  from 
all  atomic  interactions  can  enter  the  scattering  matrix 
elements  [Eq.  (2)],  and  longer-range  terms  may  be  impor¬ 
tant  in  polar  semiconductors  with  Coulomb  interactions. 
To  test  the  effect  of  including  only  the  nearest  neighbors 
described  above,  we  calculated  kl  for  Si,  Ge,  GaAs,  and 
GaP  by  using  anharmonic  IFCs  obtained  from  two 
approaches:  (i)  the  real-space  approach  described  above 
and  (ii)  a  reciprocal-space  density  functional  perturbation 
theory  approach  [19,20,29]  which  included  short-ranged 
interactions  out  to  seventh  nearest  neighbors  and  long- 
range  Coulomb  interactions.  We  find  less  than  3%  differ¬ 
ence  in  the  corresponding  room  temperature  kl  values 
from  the  two  approaches. 

Results. — Figure  2  shows  kl(T)  of  wurtzite  GaN  be¬ 
tween  100  and  500  K.  The  circles  and  triangles  correspond 
to  the  measured  kl  of  Refs.  [6,7],  respectively.  All  curves 
were  determined  by  using  the  set  2  lattice  constants,  which 
gave  the  best  fit  to  the  experimental  data  in  Fig.  1.  The 
black  curves  are  the  calculated  kl  for  isotopically  pure 
GaN,  /rpure,  in  which  all  Ga  atoms  are  69Ga.  The  red  curves 
are  the  calculated  kl  for  naturally  occurring  Ga  isotope 
concentrations,  /cnatural.  The  solid  curves  correspond  to  Km 
(in-plane  transport),  and  the  dashed  curves  correspond  to 
/cout  (transport  along  the  c  axis). 

The  phonon-phonon  scattering  rates  increase  with  tem¬ 
perature,  while  the  phonon-isotope  scattering  rates  are 
temperature  independent  [see  Eq.  (3)].  Thus  starting  from 
500  K,  with  decreasing  T,  phonon-phonon  scattering 
becomes  weaker,  causing  kl  to  rise  and  making  isotope 


temperature  (K) 


FIG.  2  (color  online).  Calculated  thermal  conductivity,  vs  tem¬ 
perature  for  naturally  occurring  (Knaturai,  two  lower  red  lines)  and 
isotopically  pure  (K-pure,  two  upper  black  lines)  wurtzite  GaN 
given  by  set  2  lattice  constants  (described  in  the  text).  Solid  lines 
correspond  to  Km  and  dashed  lines  to  Kout.  The  symbols  corre¬ 
spond  to  experimental  data  from  Ref.  [6]  (circles)  and  Ref.  [7] 
(triangles). 

impurity  scattering  more  important.  This  is  evident  from 
the  increasing  percent  isotope  effect  P  and  seen  in  the 
separation  between  the  black  and  red  curves. 

The  calculated  curves  for  wnatural  are  in  good  agreement 
with  the  experimental  data  [6,7]  and  have  room  temperature 
values  Km  =  242  W  m~ 1  K_  1  and  Kout  =  239  W  m~ 1  K_  1 . 
We  have  estimated  the  effects  of  other  point  defect  scatter¬ 
ing  such  as  from  N  isotopic  impurities,  Ga  vacancies,  and 
oxygen  (O)  substitutions,  representing  them  by  the  simple 
mass-difference  scattering  rates  from  Eq.  (3)  and  taking 
the  concentrations  reported  in  Ref.  [7]  (Ga  vacancies 
~1018/cm~3  and  O  impurities  ~102°/cm-3).  We  find 
that  including  N  isotopic  impurities  and  O  substitutions 
has  almost  no  effect  on  /cnatural  at  room  temperature,  and 
including  Ga  vacancies  lowers  /cnaturai  around  2.5%. 

The  calculated  curves  for  /cpure  are  well  above  the  mea¬ 
sured  data  and  represent  the  intrinsic  upper  limit  for  kl  in 
GaN.  At  room  temperature,  /cin  =  401  Wm-1  K-1  and 
/r0ut  =  385  W  m “ 1  K-1  corresponding  to  P  ~  60%-65% 
for  wurtzite  GaN.  With  decreasing  temperature,  this  iso¬ 
tope  effect  increases,  nearing  100%  at  T  —  100  K.  At  even 
lower  T,  sample  size  begins  to  play  a  role,  and  the  interplay 
between  phonon-phonon  and  phonon-isotope  scattering 
with  that  from  sample  boundaries  causes  a  peak  in  kl 
(not  shown).  We  have  included  phonon  boundary  scatter¬ 
ing  here  by  using  an  empirical  relaxation  time  [9]  r®  = 
L/\vk\,  with  L  —  1  mm  as  a  measure  of  the  sample  size 
consistent  with  experiment.  For  all  cases,  kl  peaks  around 
T  =  30  K  with  Kpure  =  22  000  W  m~ 1  K“ 1  being  almost  7 
times  larger  than  /<natul.al  =  3300  W  m-1  K-1  (P  ~  600%). 
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TABLE  I.  Calculated  mass  mismatch,  frequency  gap  between  acoustic  and  highest  optic 
phonon  branches,  LA  phonon  frequency  at  the  X  ( M )  point  for  zinc  blende  (wurtzite)  structure, 
^natural  at  T  =  300  K,  and  percent  isotope  effect  (P)  for  Ga-based  semiconductors. 


|  l  _  "*Ga  | 
W  anion 

Frequency  gap  (THz) 

coLA  (THz) 

^natural  (W  HI  K  ) 

P 

GaAs 

0.07 

0 

6.70 

52 

4 

GaSb 

0.43 

0.75 

4.83 

42 

8 

GaP 

1.25 

2.48 

7.70 

123 

14 

Cubic  GaN 

3.98 

6.01 

10.6 

215 

68 

Wurtzite  GaN 

3.98 

6.11 

9.34 

242 

66 

To  understand  the  above  findings  better,  we  have  calcu¬ 
lated  the  kl  for  zinc  blende  GaY  compounds  with  X  =  N, 
P,  As,  and  Sb  spanning  a  wide  range  of  properties,  given 
in  Table  I.  The  mass  variance  parameter  for  antimony 
(57.21%  121Sb  and  42.79%  123Sb)  is  gSb  =  6.59  X  10“5, 
while  the  negligible  isotope  mixtures  for  phosphorus  (P) 
and  arsenic  (As)  give  gP  =  gAs  =  0.  Figure  3  shows  P  vs  T 
for  wurtzite  GaN  and  zinc  blende  GaN,  GaAs,  GaSb,  and 
GaP.  The  smallest  P  values  occur  for  GaAs,  with  P  =  4.4% 
at  T  =  300  K.  Our  calculated  /cnatura i  and  P  for  GaAs  are  in 
good  agreement  with  kl  measurements  made  for  isotopi- 
cally  enriched  GaAs  crystals  (99.4%  71  Ga  and  0.6%  69Ga) 
which  found  a  P  =  5%  at  300  K  [30].  The  P  for  the  GaN 
polymorphs  lie  well  above  those  for  the  other  zinc  blende 
compounds. 

One  might  naively  assume  that  GaSb  would  have  the 
largest  enhancement  to  kl  with  isotopic  enrichment  be¬ 
cause  of  the  larger  isotope  mixture.  However,  for  a  given 
isotope  mixture,  it  is  the  interplay  between  the  phonon- 
phonon  scattering  and  the  phonon-isotope  scattering  that 
determines  P.  Comparing  successively  GaAs,  GaSb,  GaP, 
and  GaN  (Table  I),  we  find  that  the  phonon-phonon  scat- 


temperature  (K) 


FIG.  3  (color  online).  Calculated  percent  isotope  effect,  P,  vs 
temperature  for  GaAs  [blue  (lowest)  line],  GaSb  [green  (lower) 
line],  GaP  [red  (middle)  line],  wurtzite  GaN  (dashed  black  line), 
and  zinc  blende  GaN  [solid  black  (highest)  line]. 


tering  is  weakened  by  (a)  increasing  mass  difference  be¬ 
tween  Ga  and  anion  atoms,  which  produces  a  larger 
acoustic-optic  frequency  gap,  as  seen  in  Fig.  1  for  GaN, 
and  (b)  increasing  atomic  bond  stiffness  and  decreasing 
mass  of  constituent  atoms,  which  raise  the  entire  phonon 
spectrum  to  higher  frequencies.  Consider  point  (b)  first. 
Higher  overall  phonon  frequencies  reduce  the  phonon 
populations  of  all  modes  and  also  reduce  the  phase  space 
for  three-phonon  scattering  [31],  thereby  decreasing  the 
phonon-phonon  scattering  rates.  We  note  that  increasing 
bond  stiffness  also  increases  the  anharmonic  IFCs,  thereby 
enhancing  the  phonon-phonon  scattering  matrix  elements 
[Eq.  (2)],  but  the  stronger  influence  from  the  higher  phonon 
frequencies  ultimately  leads  to  reduced  scattering  rates 
with  increasing  bond  stiffness  in  materials  such  as  GaN. 
We  further  verified  the  importance  of  the  frequency  scale 
in  determining  kl  and  P  by  scaling  down  the  full  frequency 
spectrum  in  wurtzite  GaN.  The  resulting  increased  phonon- 
phonon  scattering  rates  significantly  lowered  kl  and  P. 
Reduced  acoustic  phonon  velocities  also  lowered  kl. 

Consider  now  point  (a).  Optic  phonons  play  a  critical  role 
in  determining  kl  of  GaN.  Although  optic  phonons  do  not 
significantly  contribute  directly  to  a  thermal  current,  they  do 
provide  important  scattering  channels  for  the  heat-carrying 
acoustic  phonons  as  has  been  pointed  out  recently  for 
silicon  [19],  diamond  [20],  and  lead  telluride  [32-34].  The 
extremely  large  frequency  gap  in  GaN  has  a  dramatic  effect 
on  kl  and  P,  because  it  strongly  restricts  three-phonon 
scattering  between  acoustic  and  optic  phonons  [31]  through 
the  energy  conservation  seen  in  Eq.  (2)  [35,36].  On  the  other 
hand,  the  isotopic  impurity  scattering  is  not  much  affected 
by  the  gap.  This  point  is  highlighted  by  comparison  of  GaN 
with  silicon.  Si  has  a  lighter  unit  cell  mass,  higher  acoustic 
phonon  frequency  scale  and  acoustic  velocities,  smaller 
anharmonic  IFCs,  and  similar  isotope  scattering  (for  natu¬ 
rally  occurring  Si)  compared  to  GaN.  These  features  sug¬ 
gest  that  Si  should  have  larger  kl  and  P  compared  to  GaN. 
However,  Si  has  no  acoustic-optic  frequency  gap.  Thus,  we 
find  that  the  calculated  /rpure  for  Si  at  room  temperature 
(155  Wm~'  K" 1 )  is  nearly  3  times  smaller  than  /cpure  for 
GaN  and  the  calculated  P  (7.6%)  for  Si  is  almost  9  times 
smaller  than  that  for  GaN  [37,38].  These  differences  stem 
primarily  from  the  large  frequency  gap  in  GaN  and  the  much 
weaker  acoustic-optic  phonon  scattering  that  results. 
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Summary  and  conclusions. — We  have  examined  the  lat¬ 
tice  thermal  conductivity  kl  of  GaN  and  its  dependence  on 
isotopes  by  using  a  microscopic,  first  principles  approach. 
The  calculated  kl  with  naturally  occurring  isotope  con¬ 
centration  is  in  good  agreement  with  measured  values  over 
a  broad  temperature  range,  and  the  increase  in  kl  with 
isotopic  enrichment  (~  65%  at  room  temperature)  is  com¬ 
parable  to  that  of  the  highest  observed  enhancements  in 
diamond,  graphene,  and  boron  nitride  nanotubes.  We  find 
that  the  large  kl  and  isotope  effect  in  GaN  result  from  the 
large  Ga  isotope  mixture,  the  mass  difference  between 
constituent  atoms,  the  light  N  mass,  and  strong  bond  stiff¬ 
ness.  The  resulting  large  frequency  gap  between  the  acous¬ 
tic  and  optic  phonon  branches  and  the  high  phonon 
frequency  scale  lead  to  unusually  weak  anharmonic 
phonon-phonon  scattering.  This  interpretation  is  validated 
by  comparison  with  GaAs,  GaSb,  and  GaP.  The  rigorous 
understanding  of  phonon-phonon  scattering  and  kl  pre¬ 
sented  here  for  GaN  gives  important  insights  into  the  origin 
of  thermal  conductivity  and  the  role  of  isotopes  in  other 
systems. 
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